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1 The general problem class 

We consider on the interval (— 7r,7r) the singular non-symmetric differential 
equation 

in which / is a 27r-periodic function having the following properties: 

/(x + 7r) = -/(x), /(-x) = -/(x); (1.2) 
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and also 



f{x)>0 forxG(0,7r) 



(1.3) 



In particular it follows that /(vrZ) = 0. We assume that / is continuous, and 
differentiable except possibly at a finite number of points, the points of non- 
differentiability excluding nZ. We assume that /'(O) = 2/7r and that < e < vr. 

Our interest in (II. ip is primarily motivated by [BeO'BSa] and |Da2] . where 
f{x) = (2/7r)sina;, and therefore (II. Ih takes the form 

d / . du\ .du ^ ,^ ^. 

le— smx— + i-r = -^^) (l-^j 
ax \ dx J dx 

with < e < 2. Eq. (II. 4p arises in fluid dynamics, and describes small os- 
cillations of a thin layer of fluid inside a rotating cylinder. It has been long 
conjectured that, despite the fact that (II. 4p is highly non-self-adjoint, all the 
eigenvalues of p. 41) are real. For a fixed non-zero e the conjecture has been 
supported by some numerical evidence, see e.g. [ChPej . who have also proved 
that eigenvalues accumulate at infinity along the real line. Davies |Da2| estab- 
lished that the spectrum is discrete, nonetheless the eigenfunctions do not form 
a basis. He also obtained a number of very interesting estimates. Finally, just 
recently Weir [We] proved, using Davies' estimates, that the spectrum of the 
highly non-selfadjoint problem (jl.4p is indeed purely real. 

The goal of our paper is to extend Weir's result and to prove that all the 
eigenvalues of (II. Ih are real under minimal restrictions on the coefficient f{x). 
One fundamental condition is the anti-symmetry of / with respect to x, see 
(II. 2p . and also some symmetry properties of the whole operator, see below. 
This makes (11.11) reminiscent (but not identical) to a wide class of so-called 
PT-symmetric equations: non-self-adjoint problems, which are not similar to 
self-adjoint, but which nevertheless possess purely real spectra due to some 
obvious and hidden symmetries. For examples, surveys, and interesting recent 
developments see e.g. [BeBol IDoDuTall IDoDuTa2l IShTl [Sh2l [LaTR IKrBiZn] 
and references therein. 

An additional difficulty in dealing with (II. ip (compared to standard PT- 
symmetric problems) is the presence of singularities at x = — vr, 0,7r. On the 
other hand, it is exactly these singularities which allow us to use a strategy of, 
very roughly speaking, matching asymptotic expansions and exact solutions. 

In the remainder of this section we shall describe some basic properties 
of the equation and its solutions. Section 2 discusses an illustrative explicit 
example in which / is piecewise linear, while Section 3 contains the general 
results. Apendices A and B cover details of some basic asymptotics and 
spectral theoretical properties of (II. ip . 
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By an elementary Frobenius analysis (see Appendix [A] for a rigorous sum- 
mary of the main asymptotic results) the differential equation (jl.ip is seen to 
possess, for each A G C, a unique (up to scalar multiples) solution in tt, vr), 
which we denote by 4>{x, A). This solution may be normalized by the condition 

<^(0,A) = 1. (1.5) 

Any solution linearly independent of (p will blow up like x"'^/'-^^-' as x ^ 0. 
Later we shall assume the more restrictive condition < e < tt/2, which is 
necessary to ensure that the differential equation has a unique solution in a 
weighted L^-space naturally associated with the problem. 

Eqn. (II. ip also has singular points at ±7r; however all solutions are square 
integrable at these points (in the unweighted space). In the weighted space it 
will turn out that whether these solutions are square integrable or not depends 
on how small e is. 

By the periodic eigenvalue problem associated with we mean the prob- 
lem of finding solutions which are square integrable and which satisfy the peri- 
odicity condition vr) = u{tt). In view of the foregoing comments about the 
solution (j) it follows that /i is an eigenvalue if and only if 

(^(-7r,/x) = 0(7r,/i). (1.6) 

Now consider, for any A G C, the function tl;{x,X) = (/>(— x,A). It is easy to 
see that ip solves (jl.lft with A replaced by —A; also, ip is square integrable and 

■ip{0,X) = l, hence ip{x, X) = (j){x, —X). In other words, 

(P{-x,X) = ^{x,-X). (1.7) 
In a similar way, taking complex conjugates in (jl.ip shows that 

^{x, A) = ^(x, -A) = 0(-x, A). (1.8) 

Eqn. (II. 8p gives a reflection principle for across the imaginary A-axis, on 
which (f) is real valued. 

Using the eigenvalue condition (II. 6p and the symmetries (jl.7p and (II. Sp . 
makes it clear that eigenvalues come in quadruples: (A, —A, A,— A). Also, the 
condition (jl.6p for /i to be an eigenvalue can be written in the equivalent form 

(p{7r,fi) = (p{7r,-fi) = (j){TT,-p). (1.9) 

We would like to show that all eigenvalues of the problem are real. 

Finally in this introduction, we emphasize that all of our results can be 
placed in a proper operator theoretic framework: in particular, the eigenvalues 
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to which we have casually referred are indeed eigenvalues of a closed non- 
selfadjoint operator. Let Cpg|.(— vr, vr) be the space of periodic functions 
in (— vr, vr). The differential expression L with domain Cpg|.(— vr, vr) defines a 
closable linear operator acting on the space L'^{—TT,7r), despite the fact that 
is not integrable at x = 0, zbvr. In order to see this we will show that L 
is a closed operator in the maximal domain 

Pm := {u G L'^{--n-,Tr) : Lu G L^(-7r,7r)}. 

Clearly Cpg|.(— vr, vr) C V^- In Appendix[B]we will show that Pm C H^{—tt, vr). 
This ensures that any eigenfunction of {L,C'^^^{'k,tt)) is a solution of the pe- 
riodic eigenvalue problem associated with (II. IH . and vice versa. Note however 
that in the general case (II. ip we do not claim to prove that the spectrum is 
purely discrete, though we believe this result is true - see Corollary IB. 21 

2 A piecewise linear example 

As an example we start with the case where / is given by 

{-2vr-i(a; + vr) -vr < x < -vr/2, 
2-n~^x -vr/2 < x < vr/2, (2.1) 

27r~i(vr-x) vr/2 < x < vr. 

In this case the equation can be solved explicitly in terms of Bessel functions. 

For X G [—vr/2, vr/2] we have 

(j){x, A) = r(z^ + l){iv\xy'^Ju{2^w\x) =: Co(x, A), (2.2) 
where v = vr/(2e). For x G [vr/2,vr] we can write 

</.(x, A) = A{\)Ci{x - vr, A) + B{\)C2{x - vr. A) (2.3) 

where 

Ci{z,\) = {iv\zy/'^Jy{2^/-iiy\z), (2.4) 

(2(z, A) = {ivXzy/^J^y{2^/-iuXz), (2.5) 

and where A and B must be determined to ensure that (j) is continuous and 
differentiable at vr/2: 

/ Ci(-vr/2,A) C2(-vr/2,A) \ / ^(A) \ / Co(vr/2,A) A 

V C;(-vr/2,A) C^(-V2,A) ) \ B{X) J [ C^(vr/2,A) J ' ^"'^^ 
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In view of the fact that Ci(0, A) = the condition (jl.9p becomes, upon using 

mi), 

B{fi)C2iO,l^)=B{-fi)C2{0,-fi). 

The properties of the Bessel functions show that C2(0, A) = i'^ /r{l — v) for all 
A, and so the eigenvalue condition becomes 

The solution of (I2.6jl is given by 



Ci(-7r/2,A) Co(vr/2,A) 
sr.. vC((-vrAA) qvrAA), 

^ ' ' Ci(-7r/2,A) C2(-7r/2,A) ^ ' ^ ' ' 



det 
det 



C((-vr/2,A) C2(-^/2,A) 

Further progress at this stage appears to rest on explicit calculation with 
Bessel functions. Firstly, it may be shown (either by reference to Watson's 
book on Bessel functions or by an explicitly solving the first order ODE satisfied 
by the Wronskian) that 

Ci(z,A) (2(2;, A) \ _ sm(i/7r) 



Secondly, observe that if we define 



x(x,A) = Jy{2\J iuXx), 

then 

Co{x, A) = T{iy + l)iii^Xx)-'/\{x, A) = cX-''/^x-'/\{x, A), 

Ci(x, A) = {iuXxr/\{-x,X) = CA-/V/2x(-x,A). 
An explicit calculation shows that 

C^(7r/2, A) = cA-^/2(^/2)-^/2(y (vr/2, A) - -x(^/2, A)), 

vr 

C((-^/2, A) = CA'^/2(-^/2)-/2(-x'(7r/2, A) - -x(V2,A)), 



I' 

IT 

V 

vr 

and hence 
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Since x'{x,X) = ivXj xJ^i^^/ iv\x), it follows that for some constant C 
independent of A, 

B{\) = CX^^^~''Jl{V2iuX7T)MV2iuXTr). 
Hence the eigenvalue condition becomes 

F{X) = F{-X), where F{X) = X^/^"" Jl{V2iuXTr)J^{V2iuXTr). (2.10) 

Setting z = \j2ivX'K, we observe that \J —2ivX'K = ±iz. 
Define 

Then a necessary condition for eqn. (j2.10p to be satisfied becomes 

\p{z)\ = 1. 

Theorem 2.1. For any non-real X (corresponding to arg{z) {ib7r/4, ib37r/4}j, 
the equation \p{z)\ = 1 cannot be satisfied. In fact, for —ir/A < arg(ibz) < 
7r/4, we have \p{z)\ < 1; for 7r/4 < arg(ib2;) < Stt/A, we have \p{z)\ > 1. 
Consequently all the eigenvalues of p. II) . with f{x) given by (12. ip , are real. 

Proof. For real v, the zeros of the Bessel function and of its derivative are 
all real. In the sector — 7r/4 < arg(z) < 7r/4, if z 7^ then iz is not real and 
so Jl{iz)Jy{iz) is nonzero. Consequently p{z) is well defined and analytic. 
Standard asymptotics show that for large \z\ in this sector, \p{z)\ < 1. On the 
rays bounding the sector, \p{z)\ = 1. By the Phragmen-Lindelof principle, since 
p \s a non-constant function which is analytic in a sector of angle 7r/2, and has a 
growth order less than or equal to two, p attains its maximum modulus strictly 
on the boundary rays. Consequently in the sector — 7r/4 < arg(z) < 7r/4, one 
has \p{z)\ < 1. 

The results for the other sectors follow similarly; for instance, if —3tt/A < 
arg(2;) < — 7r/4 then z = —it with — 7r/4 < arg(t) < 7r/4, and so \p{z)\ = 
l/\p{t)\. □ 



3 Results for the general case 

For the general case we return to our equation (|1.9p . Throughout this section 
we assume < e < ir/2. 

We start with an auxiliary result concerning the location of the zeros of 

(t>{7T,X). 
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Theorem 3.1. All the zeros of(j){'K, A) lie on the negative imaginary axis. 

Proof. Consider the differential equation p. II) satisfied by (/)(a;,A), restricting 
attention to the sub-interval (0, vr). Note that / does not change sign in this 
sub-interval. Using standard integrating factor techniques we can transform 
dET]) into 

where 

p'/p = f'/f + yief). (3.2) 
Simple asymptotic analysis shows that 

J x'^^^'^'^^ f{x) {x near zero); , , 

P^'^'^^X (n-xr^/i^^-^fix) (xnearvr). ^^'^^ 

Setting w{x) = p{x)/ f{x) and £ = i\/e, we can write p.lj) in Sturm-Liouville 
form as 

- {p{x)u)' = lw{x)u, X G (0,7r). (3.4) 

Observe that the weight function w, which is positive, has very different be- 
haviours at different endpoints: 



w{x) 



x'^l^'^'^^ [x near zero); 

(vr — x)~'^/(^^) {x near vr). 



By Frobenius analysis, there exist solutions of (11.11) with behaviours u(x) ~ 
j;-'^/(2£) and u(x) ~ 1 near a; = 0, and behaviours u(x) ~ 1 and u{x) ~ 
(vr — xyl^'^'^^ near x = vr. It is therefore easy to check that for < e < vr/2 
there is always precisely one solution in L,^ near the origin, one solution in L,^ 
near a; = vr. The problem (13. 4p is therefore of limit-point type at the two ends 
of the interval (0, vr). 

Since p > and > on (0,vr), p.4p is automatically selfadjoint in 
without the need for boundary conditions. The eigenvalues I, which must be 
real, will be precisely those values of i for which (/)(x,A), which is always in 
L\j near the origin, is in near a; = vr. This only happens when 0(x, A) ~ 
(vr — xyl^'^^^ , that is, when (^(vr, A) = 0. Since I = iXje, this establishes 
that the zeros of (/)(vr, A) lie on the imaginary axis in the A-plane. Since the 
eigenvalues of fl3.4p are all strictly positive, we can further say that the zeros of 
<\>{ti^X) all lie on the negative imaginary A-axis. □ 
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Remark 3.2. This result, which we have proved for < e < tt/2, is also true 
for < e < vr. In the case 7r/2 < e < vr the problem is limit circle in at 
X = IT but the condition (/){tt,\) = arises naturally by imposing u(7r) = in 
p. 41) . which is the Friedrichs boundary condition for this case. 

We now make the transformation A = iz'^ and define 

g{z)=cl){7T,-iz^). (3.5) 

The eigenvalue condition (^(vr, A) = (/>(7r, — A) thus becomes 

g{z)=g{iz). (3.6) 

In the A-plane, by Theorem 13.11 all the zeros of (/)(7r, A) are of the form A = —ir, 
r > 0; hence in the 2:-plane all the zeros of g lie on the imaginary z-axis and 
occur in pairs ±ian, n G N, < ai < a2 < • • • (observe that g{0) / 0). When 
z is not one of these zeros the function 



is well defined and analytic. The eigenvalue condition O.Gh becomes 

p{z) = 1. (3.8) 

Theorem 3.3. Suppose that the solution (/)(x, A) is, for each x G [0, vr], an ana- 
lytic function of \ with growth order at most 1/2. Then for —it/ A < arg(ibz) < 
7r/4, \p{z)\ < 1; for 71/4, < arg(ibz) < 37r/4, \p{z)\ > 1. Consequently, eigen- 
values of (II. ip can only lie on the lines arg(z) = 7r/4 (mod vr/2), which, in 
terms of X = iz"^ , means that all eigenvalues are real. 

Proof. We start by outlining the idea of the proof; after the outline, we fill in a 
technical detail. 

Firstly we know that g{^) ^ 0, that g is an analytic function of 2^ and 
that the zeros of g occur in pairs ±ian. If these zeros have the appropriate 
asymptotic behaviour then we can write 



00 y 2 \ 

eMh{z'))ll(l + ^) (3.9) 

n=l ^ 



where h is an analytic function and the infinite product is convergent. Since 
we also know that g has exponential growth order at most 1, it follows that h 
must be constant, and hence 



71=1 
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It is then a simple matter to check that for each n G N, 



1 - 


2^ 






1 + 


"1" 



< 1 Vz such that - 7r/4 < arg(z) < 7r/4, 



and the result is immediate. 

It only remains to check that the zeros of g grow sufficiently rapidly to 
ensure convergence of the infinite product in p.9p . To this end we recall from 
the proof of Theorem I3.1l that the zeros of (/)(7r, A) are given by the eigenvalues 
of the ODE (13. ip on (0, vr) with Friedrichs boundary conditions. We therefore 
wish to examine the large eigenvalues of this problem. To this end we make the 
transformation 



{x,X) 



1 



--'4'{x,x). 



where p is given by (I3.2p . This reduces the differential equation to 



(3.10) 



(3.11) 



where Q = p~^/^(p~^/^p')'/2 and where /i = ie^^X is already known to be real 
and positive. A simple calculation based on the known asymptotics of p near 
a; = and x = vr shows that 



Q{x) 



4 \ 2e 



-1) X '■ 



X 
X 



0, 
vr. 



Since < e < 7r/2, we have Q{x) —>■ +oo for x \ and x y ir. By 
Sturm Comparison we can therefore obtain lower bounds on the eigenvalues /i 
by discarding Q and considering the eigenvalues fi of the problem 

fix) 

where Q > \s a sufficiently large constant. Standard WKB estimates (see, 
e.g., Bender and Orszag [BeOr]) give for the nth eigenvalue the formula 



fix) 



dx ~ n, 



where we observe that the integral converges because / has simple zeros at the 
endpoints of the interval. Hence 



flri = 0{7?). 
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Recalling the transformations A = iz^ and A = ie~^fi and the definition of the 
zeros On of g and their correspondence with the eigenvalues of this problem, 
we must have 

Aj2 is ^72 5 

whence 

an = 0{n). 

This is sufficient to ensure convergence of the product in O.Qp . □ 



Remark 3.4. In the case of the piecewise linear coefficient example of eqn. (12. ip . 
(j){x,\) has growth order 1/2 in A and hence 0(x, —iz^) has growth order 1 in 
z, for any x e [0,7r]. This follows from the explicit expressions for the solutions 
given in terms of Bessel functions in Section [2j 

The same conclusion also holds for the case f{x) = (2/7r) sm(2;) considered 
by Davies [Da2j : the equation admits a change of variables t = tan(x/2), and 
the solutions are then given in this case in terms of Heun G-functions on the 
interval [0, +oo). For relevant asymptotic results, see e.g. [Fi] . 

In both of these cases, therefore, the eigenvalues of the corresponding PT- 
symmetric eigenvalue problem are all real. 

The following lemma widens the class of functions / for which the function 

(f){x,\) has growth order 1/2 in A for all x G [0, vr]. 

Lemma 3.5. Suppose that f falls within the class introduced in Section{^and 
that f is linear in arbitrarily small neighbourhoods of x = 0, x = n. Then 
4>{x, A) has growth order 1/2 in A for all x G [0, vr]. 

Proof. Suppose that f{x) = 2x/it for x £ [0, 5] and that f{x) = 2{'k — x)/-k for 
X G [vr — (5, vr] for some (5 > 0. Then the explicit expression for (j) given in ()2.2p 
is valid for all x € [0, (^] and, with standard asymptotics of Bessel functions, 
establishes that (f){x, A) and (l)'{x, A) have exponential growth order 1/2 in A for 
all such X. Writing the differential equation for 4> as 

±( 4>{x,\) \^( Ai/V/(x)\/ <^(x,A) 
dx V X~^'^f{xW{x, \)J V -iX'^Ve -l/(e/) J { X~y^f{x)^'{x, A) 

(3.12) 

we observe that this is a system of the form Y' = A{x, X)Y in which \\A{x, A)|| < 
C^jAj^/^ for all x £ [5,tt — 6]. By standard Picard estimates 



0(x, A) 
X-y^f{x)ct>'{x,X) 



< 



<P{6,X) 
X-y^f{5)ct>'{S,X) 



exp(Q(x-5)|A|i/2) 
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for all X G [6,7r — S]. This establishes that the result of our lemma holds for all 
xG[0, vr — 5]. Finally we must extend the result over [vr — 5, vr]. To this end we 
write 

4>{x,X) =A{X)Ci{x-Tr,X) + B{X)C2{x-Tr,X), x e [tt - 6,7r] (3.13) 

where Ci and ^2 are the functions introduced in eqns. (j2.4p . (j2.5l) , both of 
which have growth order 1/2 in A for all x. The coefficients may be found by 
solving the system 

/ CiM,A) C2{-S,X) \fA{X)\_f </.(<5,A) \ 

V p{5)C[{-6,X) p{5)C^{-6,X) J [ B{X) J - [ p{6)4>'{5,X) J ^^■'^> 

where p may be any function which is non-zero at 6 but is most conveniently 
chosen according to p.2jl . This ensures that the determinant of the matrix on 
the left hand side of (j3.14p is independent of 6 and may be calculated from its 
value as 6 —>■ 0. From the asymptotics of Ci and C2 it turns out that this value 
is a non-zero constant, and so 

( AW A _ ( P{S)C'2{-S,X) -C2{-S,X) \ ( <P{5,X) \ 

\ B{X) ) ^^^^ V -P('5)C((--5,A) Ci(-5,A) ) V Pm'{5A) ) 

Every quantity appearing on the right hand side has growth order at most 1/2 
in A and so ^(A) and B{X) have growth order 1/2 in A. Thus from p.l3p it 
follows that 4){x, A) has growth order 1/2 for all x G [vr — 5, vr]. □ 

Remark 3.6. The ideal way to obtain our final result now would be to show 
that (j){x,X) always has growth order 1/2 for the whole class of coefficients / 
introduced in Section [H Unfortunately this does not appear to be easy, since 
growth orders are not necessarily preserved in approximation limits. For instance, 
exp(A^) has growth order 2 but can be approximated locally uniformly by Taylor 
polynomials, which all have growth order 0. Nevertheless the following result is 
true. 

Theorem 3.7. Let f be a function of the class introduced in Sect/on [I] wliich 
additionally has the property that f"{0) and /"{it) exist. Then all the eigen- 
values of eigenvalue problem p.ip are real. 

Proof. By Theorem 13.31 and Lemma [33] the result holds for functions / which 
are linear in arbitrarily small neighbourhoods of the endpoints. 

Now suppose that we have a problem for which / is not linear near the 
endpoints. Nevertheless we may approximate / by a function which is linear 
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on [0,6] and [vr — 6,tt] for small 6. From (11.60 and (jl.7p , we know that the 
condition for /i to be an eigenvalue is that /i be a zero of the function 

d{X)=(l){Tr,X)-(l){Tr,-X). 

We shall attach subscripts S to the quantities where / is replaced by a linear 
function on [0, 5] U [tt — 5, vr] : thus we shall approximate A) by (j)s{x, X) and 
d{X) by d5{X). If we can show that 

VunidsiX) - d{X)) = (3.15) 

locally uniformly in A, then the zeros of d{X) will be precisely the limits of the 
zeros of ^^(A). This is achieved, for instance, by using the argument principle to 
count zeros inside any contour, exploiting the fact that for analytic functions the 
locally uniform convergence of the functions implies locally uniform convergence 
of their derivatives. Since the zeros of ds all lie on the real axis, we shall have 
proved our result. 

Choose fs as follows. On [0, (5] U [vr — 5, vr] fs should be linear and it should 
also match the values of f at 0, 6, ir — 6 and vr: consequently we shall have 

/5(0) = /(0)=0; fm=fi6); /^(vr - 5) = /(vr - 5); /^(vr) = /(vr) = 0; 

also, as (5 \ 0, 

- no); /Kvr) - fin). 

Frobenius analysis (more precisely, the uniform asymptotics in Lemma IA.4[1 
shows that 

^{x,X) = l- ^^fj^^^^^ x + oix); ,^'(x,A) = -^-p^^ + o(l) (3.16) 

for small x. It is convenient also to have the second solution <I>(x,A) which 
satisfies 

$(x, A) = x-^/("-^'W)(l + 0(x)), $'(x, A) = ^^x-i-i/(^^'W)(l + 0{x)) 

(3.17) 

for small x. 

Note that we also have similar asymptotics for (ps and $5: 

Mx, A) = 1 - YTJflio)'' + ^) = - i + fm + 

^>(x,A) =x-i/("-^^(°»(l+0(x)), $'(x,A) = -^x-i-i/("^«(°))(l+0(x)). 
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Moreover the correction terms can be bounded locally uniformly in x G 

by the remarks following Lemma lAA\ Consequently, since /^(O) — /'(O) = o(l), 

cP{6, A) - MS, A) = 0(6); cl>'i6, A) - <p'siS, A) = o(l). 

All of these bounds are locally uniform in A. 

Now on [6, TT — 6] we have 

(J)s{x,X) = ci(l){x,X) +C2Hx,X) (3.18) 

where the coefficients ci and C2 are chosen to ensure continuity and differentia- 
bility atx = 6. In terms of the coefficient p introduced in 0.2jl . and suppressing 
the A-dependence for simpler notation, we have 

/ m m \( = ( MS) \ 

V pcP'{5) p«J>'(5) J\c2 J V pMS) J ' 

whence, in terms of the determinant A of the matrix on the left hand side of 
this system, 

1 \ / <^'is){Ms) - m) - ms) - <p'is))m \ 
; + A V mms) - cp'{5)) - ^'{6){ms) - m) ) ■ 

Bearing in mind that from (13.31) . 

and that A is a (^-independent constant, we obtain by elementary estimates that 

C, = 1+0(6); C2 = o(dl + l/(^^'(0))). 

From p.lSp it follows that for each fixed x £ [0,tt — 6] we have 

lim (j)s(x, A) = (j)(x, A) 
<5\o 

locally uniformly in A. We need to extend this convergence to x = tt. 

To this end we introduce the solutions ■ip{x) and ^(x) for the unperturbed 
system and ips{x) and ^^(x) for the perturbed system, determined by the 
asymptotic behaviours (see Lemma |A3] again) 

^(x) = 1+- 'A^{n-x)+o{Tr-x), ^{x) = (^-x)i/(^^'W)(l+0(7r-x)), 

(3.19) 
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^5^^) = 1+ 1 '),,^A 7r-x)+oi7r-x), ^s{x) = (7r-x)i/(^^^W)(l+0(7r-x)). 

(3.20) 

We write 

(psix) = ci^si'six) + C2,5'ifs{x) 

and 

= ClV'(x) + C2^'(x). 

In view of the asymptotics it follows that 

0(7r) = ci; (^^(vr) = ci,^. 

Therefore in order to establish that lim^^o 'Psi'^^ ^) = ^) locally uniformly 

in A it is sufficient for us to prove that ci 5 converges locally uniformly to ci. 
Now 

pVK^ - S) p^'sin -6) J\ C2,5 ) V - ^) 

and a similar equation holds with (5 = 0. A calculation similar to the ones 
performed before shows that 



Cl,5 = Ci + |-{(^^-^')05 + *'( 



A 



We now estimate the various terms, bearing in mind that from O.Sp we have 

p(vr-5) = 0(5i-WW)) 
which may be large. Consider, for example, the term 

p{tt - 5){%{tt -5)- ^'{tt - 6))(l)siTT - 6). 
In this term we have 

_ 5) _ vl;'(^ -S)= 0(5l/(e/'W)-l(<5l/{e/^-W)-l/(e/'W) _ 1)). 

The multiplicative factor of (5^/('^-^'('^))^^ cancels with the p(7r — 5) term, while 
the (I)s{tt — 6) term is close to (/>(7r — 6) which remains bounded as 5 \ 0, since 
all solutions of the unperturbed equation are bounded in a neighbourhood of tt. 
This leaves a term 

0{6^ - 1) 
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in which 

However the fact that / is twice differentiable at vr ensures from the construction 
of fs that 

n{7r)-f\n) = p"{7T)+oil) 
and so the term which is 0(6'^ — 1) is actually 

0(55/"W/2 _ 1). 

This tends to zero as S tends to zero, regardless of whether /"(vr) be positive 
or negative. 

The other terms may be dealt with in a similar way. Since all the asymptotics 
are locally uniform in A, we get 

lim(/>5(7r, A) = (I){it,\) 
locally uniformly with respect to A. The result follows. □ 

A Appendix: Asymptotics 

In this Appendix we examine the asymptotic behaviour of solutions u of p. II) 
and their derivatives, separately near x = and x = n. In order to have a 
completely rigorous treatment in each case we shall transform the problem to 
an infinite interval and invoke the Levinson Theorem [Eal Theorem 1.3.1]. 

Consider first the behaviour in a neighbourhood of 0. Make the change of 
variable x = exp(— i) so that maps to oo. Eqn. (jl.ip now becomes 

ie—[e^f{e~^)—]-i— = \e~\. A.l) 
dt\-'^'dt)dt ^ ' 

Under our hypotheses on / we know that for all a; in a neighbourhood of 0, 

/(x) =x/'(0)+r(x), 
where |r(x)| < Cx^ for some constant C > 0. As a consequence 

eV(e-*) = /'(0) + p(t), 

where 

|p(t)| <Cexp(-t). (A.2) 
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Consequently we can write (jA.ip as a first order system 



d 



u \ _ ( " /'(0)+p(t) 

e ^ f'{0)+p{t) 



dt \ {f'{0)+p{t))u' J \ ^e-* ^rr^ I V (/'(0)+p(t)K 



(A.3) 
or as 

^ = {A^+R{t))n, (A.4) 

where 

u 



u 



Ac 

and 



(/'(0)+p(t))n' 
l//'(0) 

e-V/'(o) 



\\R{t)\\ < Ce'K (A.5) 
The matrix A^o has eigenvalues and (e/'(0))~^ and is diagonalizable: 

A^ = VAV~\ 

where A = diag(0,(e/'(0))~i), 



V 



1 e 
1 

Consequently the system (]A.4[) can be diagonalized by the transformation v = 
V^^u, and 

^ = {A + V-'R{t)V)^r. (A.6) 
Lemma A.l. Suppose that T is chosen sufficiently large to ensure that 



\\V-^R{t)V\\dt < 1/2. 



T 



Then the system ( f/4.6j) has solutions vi and V2 stvc/j t/?at 



vi(t) = Q+ri(t), (A.7) 
V2(t) = exp(f/(e/'(0))) I (J) + r2(i)} , (A.8) 
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where, for all t > T, 

\\VR{s)V~^\\ds. (A.9) 

Consequently, in view of ( f/4.5]) and the transformation v = V~^u, the linear 
system I^A.4\) possesses solutions ui and U2 such that 

ui(t) = i^Q +m{t), (A.io) 

U2(i) = exp(t/(e/'(0))) +r?2(t)} , (A.ll) 

where, for some constant C > and for all t > T, 

max(||r?i(t)||, Hr^sWH) < C7exp(-i). (A.12) 
Proof. This is just the statement of Levinson's Theorem for our problem. □ 

It is now a simple matter to translate this lemma into information about 

u{x), u'{x) and f{x)u'{x) for x in a neighbourhood of 0: 

Lemma A. 2. For the equation there exist solutions ui and U2 such that 

ui{x) \ _ ( ^ 



f(\n\ ] - \ n ; + ^i(^)' (A.13) 
-j{x)u\{x) J \ J 

where ||?7i(a;)|| < Cx, \\rj2{x)\\ < Cx for all < x < xq, where xq is some 
positive real number. Moreover the error bounds are locally uniform in X and 
depend only on X, e and supo<a;.<^,j [''(a^)!- 

Observe that the solution ui is the solution cj) which is regular near zero, 
used throughout our article. 

A similar lemma can be proved concerning behaviour near x = tt: 

Lemma A. 3. For the equation there exist solutions and such that 



where \\r]iix)\\ < C(7r — x), \\r]2ix)\\ < C{'k — x) for all tt — xq < x < tt, where 
Xq is some positive real number. Moreover the error bounds are locally uniform 
in X and depend only on X, e and sup„^p^^^o 
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These results are sufficient for most of our article, but not quite for the 
estimates in the proof of Theorem 13.71 In particular, although Lemma lA^ gives 
the results p.l7p and p.20p . it does not quite give the results p.l6p and p.lQp . 
since the terms i:f^jr^x are absorbed into the error terms. The solution to 
this problem relies on using an integral equation, eqn. (1.4.13) in [Eal p. 12], 
according to which the solution vi of Lemma [ATI] will satisfy 

vi(t) = ( J ) - Fit)F~\s)V~'Ris)Vviis)ds, (A.17) 

in which, for our system, 

^(«)=(J exp(t/("r,0))))' '■^''^ 

The higher quality approximation is achieved by replacing vi in the integral on 
the right hand side by its approximation from Lemma IATI] so that 

V-^R{s)V^rl{s) = V-^R{s)V ( J ) + C'(exp(-2s)). (A.19) 

Suppose that /"(O) and /"(vr) exist. Some explicit calculations show that 

''-'^wv ( _.._r,o;W(o)0-io ;,)) 

where 

g{t) = 0(e*(eV(e-*) - /'(O)) - f"{0)/2) (A.20) 
tends to zero as i — > oo. Hence 

V''R{s)VMs) = e-'^Xi ) + 0{\g{s)\e-'^) + 0{e-^'). (A.21) 

Substituting (jA.lSD and (jA.2ip back into (jA.17P and transforming back to the 
original variables yields the following result: 

Lemma A. 4. The solution ui of (fJ. J|) mentioned in Lemma lA^ has an asymp- 
totic form 

where \\mix)\\ < C{x'^ + \x-'^fix)-f'{0)-xf"{0)/2\) for all x in some interval 
[0, xq], xq > 0, and the constant C depends only on A, e and supo<a,<^.g 
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As an immediate corollary of this Lemma, and in particular of the form of 
the error term, if one changes / in a neighbourhood [0,6] to some function fs, 
then the same asymptotics will hold in [0,6] with /^(O) replacing /'(O) on the 
right hand side of (jA.22p . A similar reasoning works near x = n and gives a 
rigorous justification of the asymptotic results used in the proof of Theorem [T71 



B Appendix: Linear operator realisation of L 

The quadratic form associated to L with domain C^(— 7r,7r) is strongly indefi- 
nite. In fact, its numerical range is the whole of C. In this appendix we show 
that, despite of this strong indefinitness, L is a closable operator on C^(— vr, vr). 
The key to this result depends on an analysis of the behaviour of n G Vm in 
neighbourhoods of 0, zbvr. 

Lemma B.l. Let < s < tt/2. There exists a constant Cg > 0, such that 

/TT piT 
\u'{x)\'^dx<ce \Lu{x)\'^dx (B.l) 
-TT ■/ —TV 

for all u G V^. 
Proof. Suppose that 

- iLu = e{fu)' + u' = ve L^{-^, vr). (B.2) 
We will only show that 

"TT" /'TT" 

u' {x)\'^ dx < Cs I \v{x)\^ dx. 

JO 

If we restrict our attention to x € [0, vr], the left side of (IB.2P can be written as 

[pu ) = —V, 

where p is given by (13.21) and satisfies the asymptotes (13. 3 j) . Note that p{x) 
as well as f{x) are continuous and positive in (0, vr). By virtue of Hardy's 
inequality, we can find a constant > such that 

|2 



"^'l U M2^ ^ \p{x)u'ix)\ , 

\u'ix)\'dx<a, LA_LA_l^dx 



- Jo x^/^ 

r'^ \p{x)/f{xr\v{^r . r'\(^\-A 

< ag ; dx = Ug / fix dx. 

~ Jo x^/^ Vo 



19 



Similarly, if we change variables y = tt — x, we can find be > such that 

Jo Jo y^-""/^ Jo 

This ensures the desired inequality. □ 
In the proof of this lemma observe that Hardy's inequality 

t e-^\w{t)\^ dt < f\''\w'{t)\^dt 
4 Jo Jo 

holds for any a G M if u; is absolutely continuous and supp(t(;) C (0,6). See 
[Pali Section 5.3]. 

The fact that L is closed in the maximal domain is a straightforward 
consequence of this lemma. 

Corollary B.2. Let L denote the operator defined by taking the closure of the 
differential expression at the left side of (II. ip in Cper(— vr, vr). If the resolvent 
set of L is non-empty, then the spectrum of L is the set of eigenvalues of the 
periodic problem associated with (jl.ip . 

Proof. Suppose that the resolvent set of L is non-empty, and let A be a member 
of this set. Then, by Lemma IB. II (L — A)~^ is a bounded operator mapping 
L^(— TT, vr) into H^{—tt, it). As H^{—ir, tt) is compactly contained in vr, vr), 
then necessarily (L — A)^^ is a compact operator. The spectrum of L therefore 
consists of isolated eigenvalues of finite multiplicity. These are precisely the set 
of eigenvalues of the periodic problem associated with (jl.ip . □ 

Note that if f{x) = {2/it) sin(x) the hypotheses of this corollary are satis- 
fied. See [Da2] and [Wi] . 
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